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We measured the THz reflectance properties of a high quality epitaxial thin film of the Fe-based 
superconductor BaFei.84Coo.i6As2 with Tc=22.5 K. The film was grown by pulsed leiser deposi- 
tion on a DyScOa substrate with an epitaxial SrTiOa intermediate layer. The measured Rs/Rn 
spectrum, i.e. the reflectivity ratio between the superconducting and normal state reflectance, sug- 
gests the presence of a superconducting gap A a close to 15 cm~^. A detailed data analysis shows 
that a two-band, two-gap model is necessary to obtain a good description of the measured Rs/Rn 
spectrum. The low-energy Aa gap results to be well determined (Aa=15.5±0.5 cm~^), while the 
value of the high-energy gap As is more uncertain (As=55±7 cm~^). Our results provide evidence 
of two electronic contributions to the system conductivity with the presence of two optical gaps 
corresponding to 2A/kTc values close to 2 and 7. 



PACS numbers: 74.78.-w,74.25.Gz,78.30.-j 

I. INTRODUCTION 

Since the early days of the discovery of superconduc- 
tivity in the iron-pnictide family [ij, analogies and dif- 
ferences with copper-oxide high-Tc superconductors have 
been extensively discussed. In particular, the importance 
of electron correlation Q and the degree of itinerancy of 
the charge carriers 0] is still a matter of debate. Closely 
connected to that issue, is the problem of the symmetry 
of the superconducting order parameter. Indeed nodal 
gap symmetries are more liable to accomodate Cooper 
pairs in the presence of strong repulsions, than a fully 
gapped s-wave symmetry. Unfortunately, the determi- 
nation of the gap properties of pnictides is further com- 
plicated by the fact that multiple bands are crossing the 
Fermi level, thus forming semi-metallic electron and hole 
pockets. While both nodal 0-01 and node-less [j-fTsI 
gaps have been proposed in literature, there is now a 
growing experimental evidence favouring the s-wave sce- 
nario. This finding is consistent with unconventional su- 
perconductivity mediated by spin-fluctuations , through 
the so-called s± pairing state [1J|. In this scenario, it is 
of basic importance a detailed knowledge of symmetry, 
number and values of the possible superconducting gaps. 

In studying the superconducting Fe-based compounds, 
the 122 family obtained by doping of the parent com- 
pound AFe2As2 (A=Ba, Sr, Ca) is particularly suitable 
to explore the nature of the superconducting gaps, as it 
can be prepared in sizeable single crystals [15l| as well as 
in high quality epitaxial thin films [l6l . Il7l|. Moreover, 



compounds of the 122 phase are chemically and struc- 
turally simpler and more isotropic than those of the 1111 
family (ReFeAsO, Re=rare earth), and can have a fairly 
high transition temperature. In particular, Tc can be as 
high as 25K in the electron doped BaFe2-a;Co2:As2 com- 
pound 18], which is the object of the present work. The 
superconducting gaps in this system have been measured 
by various techniques, namely angle resolved photoemis- 
sion spectroscopy (ARPES) M , heat capacity [iJI, scan- 
ning tunneling spectroscopy [IQI [20|. penetration depth 
probed by ^SR |10], point contact spectroscopy |9], Ra- 
man scattering jj], and thermal conductivity [13] ■ The 
results are widely scattered and even in sharp contrast 
with one another in terms of symmetry of the gaps, 
their number and their values. Indeed, the gaps have 
been found to be either nodal [3-0| or nodeless [sl-fisj. 
isotropic [1, Q or anisotropic [3.Tl3|. As to number and 
values of the gap(s), both a single gap [i,|l9, 20] and two 
gaps d, [13, El have been reported in the range from 
2A/KTc ~ 2 to ~ 7. It is important to note that caution 
must be used in comparing results obtained on different 
samples since doping [s^ and disorder [l^ [2ll [2^ have 
been indicated to be likely responsible for the apparent 
disagreement of literature data on gap values. Caution 
must be used also in comparing results obtained through 
different techniques since impurity effects, sample inho- 
mogeneities, and surface off-stoichiometry can differently 
affect bulk and surface sensitive measurements. 

In general, infrared spectroscopy is a powerful tool 
to study the properties of a conducting system since 
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the optical response of the free-charge carriers is deter- 
mined by the frequency-dependent complex conductivity 
(T = (Ti + ia2- In a simple one-band system, the stan- 
dard Drude model (with parameters plasma frequency 
and scattering rate 7) describes the complex conductiv- 
ity CTAT in the normal (N) state 23]. Here it is worth to 
recall that the d.c. conductivity aQ—(Ji{uj = 0) is pro- 
portional to ^2^/7. In the superconducting (S) state, 
the standard BCS model (Mattis-Bardeen equations [13] , 
with parameters ctq ^nd superconducting gap A) can de- 
scribe the complex conductivity as [H]. On this basis, 
far-infrared/terahertz measurements can be of particular 
importance since a mark of the superconducting gap A 
can be observed at huj ^ 2 A (optical gap) for an isotropic 
s-wave BCS superconductor. For a bulk sample, in par- 
ticular, a maximum at the optical gap is expected in 
the ratio Rs/Rn, where Rs and Rn are the frequency- 
dependent reflectances in the superconducting and nor- 
mal state, respectively 

Infrared spectroscopy has been widely employed in 
studying Fe-based compounds [l^ll^]. In particular, in- 
frared/terahertz measurements on BaFe2-a;Coj;As2 sam- 
ples |28l - [34j have been recently employed to study the fre- 
quency dependent conductivity and to get information 
on the superconducting gap(s). However, there is still 
poor consensus on the value and, most importantly, on 
the number of gaps resulting from the analysis of THz 
data, as well as from the analysis of data from differ- 
ent techniques, as discussed above. In the present work 
we investigated the multiple-gap character of the Co- 
doped BaFe2As2 system by performing reflectance mea- 
surements at terahertz (THz) frequencies on a high qual- 
ity film [1^ . The results of the previous THz studies will 
be reported in the following, together with the results of 
the present work. 



II. EXPERIMENTAL RESULTS 

Epitaxial Co-doped BaFe2As2 films with high transi- 
tion temperature were recently obtained by using novel 
template engineering [16]. The investigated 350 nm thick 
BaFei.84Coo.i6As2 (BFCA) film was grown in-situ by us- 
ing pulsed laser deposition (PLD) on an (110) oriented 
DyScOs (DSO) single crystal substrate with an interme- 
diate layer (20 unit ceU) of SrTiOs (STO). Out-of plane 
9-26 X-ray diffraction was used to evaluate the quality 
and the orientation of the films: The full-width at half- 
maximum (FWHM) of the 004 rocking curve is indeed 
narrow (0.55°) and the 001 peaks, visible in the spectra, 
indicate that the film is grown with the c-axis normal 
to the substrate. For further informations on the film 
growth and on their epitaxial quality see Ref.[16]. 

The temperature-dependent resistivity p{T) of the in- 
vestigated film is reported in Figll] showing the onset 
of the superconducting transition at Tc= 22.5 K with a 
width of around 2 K. The error bar in the value of the 
residual resistivity po (0.15±0.01mJ7cm) is mainly due to 
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FIG. 1: (Color online) Resistivity p{T) of the 
BaFeo.84Coo.i6 As2 film. In the inset, a zoom of the su- 
perconducting transition. The residual resistivity is po=p(25 
K)=0.15±0.01 mficm. 



uncertainty in the film thickness t estimate (350±25 nm). 

We measured the R{T)/Rm spectrum (where Rm = 
R{25 K)) of the BFCA film in the THz region by em- 
ploying synchrotron radiation at the infrared beamline 
SISSI [Hi] at the synchrotron Elettra (Trieste, Italy). We 
remark that the large size of the film surface, together 
with the use of a high-fiux synchrotron source, enables 
addressing the reflectivity at THz frequencies with supe- 
rior signal to noise ratio with respect to measurements 
on single crystals. Indeed, our previous Rs/Rn mea- 
surements on an ultra-clean MgB2 film, performed in the 
THz region by using synchrotron radiation, provided the 
first evidence of the effect of the two gaps in optical mea- 
surements [3^, not observed in single crystals 3^. The 
R{T) I Rm measurements were made by cycling the tem- 
perature in the 6-25 K range, without collecting refer- 
ence spectra, in order to avoid any variation in the sam- 
ple position and orientation, which may yield frequency- 
dependent systematic errors in R(uj). The obtained re- 
sults are reported in Fig. [5^. We remark that the high 
accuracy of the present synchrotron measurements allows 
for the detection of small effects of the order of 0.1 % at 
THz frequencies. 

For T ~ Tc, the R{T)/Rn curve is fiat within ex- 
perimental uncertainties, while on decreasing tempera- 
ture the R(T) / Rn spectrum increases on decreasing fre- 
quency until the Rs/Rn (with i?s=R(6 K)) reaches a 
plateau at about 30 cm^^. At even lower frequencies, 
Rs/Rn becomes nearly constant. The frequency depen- 
dence of Rs/Rn first suggests the presence of a supercon- 
ducting gap A close to 15 cm~^. Indeed, for w — >■ 0, the 
reflectance Rn of a conducting system tends to 1 for a 
bulk system, and to a slightly lower value for a thin film, 
since the penetration depth becomes larger than the film 
thickness. Therefore, since Rs approaches 1 at w = 2A, 
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FIG. 2: (Color online) (a) R{T)/Rn spectra of the BFCA film 
in the THz region, (b) The 6, 12, 18 K spectra conveniently 
shifted are compared with the two-band (full line) best fit 
profiles. In the 6 K case, also the one-band best fit profile 
(dotted line) is shown. In the inset of (a), the R/Rn best-fit 
model spectrum of the film is compared with the correspond- 
ing spectrum of a bulk sample (see text). In the inset of 
(b), the normalized A(r)/Ao values are reported {A a, full 
squares, As, open circles), with Ao=A(6 K). The full line is 
the standard BCS behaviour. 

Rs/Rn exhibits a maximum around 2 A in the case of a 
bulk sample, while it remains nearly constant below 2A 
in the film case [36]. 

III. ANALYSIS AND DISCUSSION 

For a detailed analysis of the measured spectra, we 
used the procedure successfully employed in the MgB2 
[3^ and VsSi cases, o-jv is described by the Drude 
model, and as by the Zimmermann model (39j which 
generalizes the standard BCS model [13, to arbitrary 
temperatures and to systems of arbitrary purity (i.e., to 
arbitrary 7 values), as thus depends on the reduced tem- 
perature t = T/Tc and on the parameters Q, 7, and A. 
When two or more bands are assumed to contribute to 
the film conductivity, we employed the parallel conduc- 
tivity model (in the two-band case jUjHj], = aA+a's), 
in both the normal and superconducting states. 

From the calculated a, by using standard relations [23| . 



it is then possible to compute the model refractive index 
n = n + ik oi the film. The reflectance spectrum of the 
three-layer system (film-STO-DSO) can then be evalu- 
ated by using the thickness and the n, k values of each 
layer [401. For STO, detailed n and k values in the entire 
infrared region are reported in the literature [4l[. Since 
such data are not available for DSO, we performed abso- 
lute reflectance measurements on a bare DSO substrate 
over a broad range, thus obtaining the n and k values 
relevant to this work. A detailed study of the DSO opti- 
cal properties in the far-infrared region will be reported 
in a forthcoming paper. We verified that, for any reason- 
able choice of the Drude parameters determining n and 
k of the conducting film, the intermediate STO layer has 
no appreciable effect on the THz refiectance spectrum 
because of its extremely small thickness. 

When only one band is supposed to contribute to the 
film conductivity (one-band, 1-b model), the unknown 
parameters which determine the Rs/Rn spectrum are 
il, 7, and A. We thus employed a three-parameters fit- 
ting procedure, and repeated the procedure for different 
values of the film thickness (<=350±25 nm). This proce- 
dure does not provide an at least satisfactory description 
of the Rs/Rn data independently of the t value, as shown 
by the best-fit curve in Fig. [^b, obtained for i=350 nm 
with n = 3.2 eV, 7=800 cm-\ and A=32 cm-\ The 
failure of the best fit profile in reproducing the experi- 
mental spectral shape, both above and below 2A, indi- 
cates the inadequacy of the 1-b model in describing the 
optical response of the film at THz frequencies. 

In a two-band scenario, the six parameter ($7^,7^, A^, 
i = A, B) fit provides a very good description of the 
Rs/Rn spectrum, as shown in Fig. ^jp. We remark 
that only a two-band model can describe the measured 
Rs/Rn spectrum, and in particular the smoothly de- 
creasing slope above 2A. The fit provides, for i=350 nm, 
plasma frequencies il-A — 1 eV and 0^ — 2 eV, — 200 
cm"-'^, 7b — 2500 cm~^, A^ ~ 15 cm~^, and Ab — 55 
cm^^. The Rs/Rn best-fit model spectrum is compared 
in the inset of Fig. O) with the model Rs/Rn spec- 
trum evaluated with the same parameters in the bulk 
case (i.e., when t — > 00). As noted above, the film spec- 
trum is nearly constant for frequencies lower than the 
peak frequency of the bulk spectrum. 

In discussing the fit results, we first remark that the 
Drude parameters f2i and 7.^ are affected by rather large 
uncertainties (of the order of 20% when the uncertainty in 
the film thickness is considered) since these parameters 
are closely inter-related. However, the impact of these 
uncertainties is reduced in the total plasma frequency 
{VLtot = \/^A + = 2-2 ± 0.3 eV) and even more 
in ao—a^A+(ToB =(7.0±0.5)10'^ f2~^cm^^, corresponding 
to a residual resistivity /cio=l/o'o=0-140±0.015 mficm, in 
good agreement with the measured value /9o=0.15±G.01 
mr2cm(see figH]). This result supports the validity of the 
employed procedure. As to the gap values, we can safely 
pose Aa=15.5±0.5 cm~^ and Ab=55±7 cm~^. Indeed 
A. A is well determined since, as noted above, Rs/Rn is 
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FIG. 3: (Color online) Contributions of A and B bands to 
o"! as evaluated by using the 2-b best fit parameter values in 
both the normal (N) and superconducting (S) state, (b) total 
(Tijv and (Tis. 



nearly constant below 2A^, while the value is much 
more uncertain since this gap does not give a clear spec- 
tral feature. Nonetheless, as discussed above, the B-band 
contribution needs to be included in the model in order 
to obtain a good description of the experimental Rs/ Rn 
spectral shape. 

We then attempted a two-band modelling of the 
R{T)/Rm spectra measured at 12 and 18 K by using 
the 6 K values of the Drude parameters and using A^i 
and As as free parameters in the fitting procedure. As 
shown in Fig. [2J5, it is not possible to obtain a satisfac- 
tory description of the measured spectra, in particular at 
low frequencies. As a consequence, the resulting /S.a and 
As values are affected by large uncertainties, as shown 
in the inset of Fig. [2b where the normalized A(T)/Ao 
values are reported, with Ao=A(6 K). In the same fig- 
ure, the standard BCS A(T)/Ao behaviour is reported 
for comparison. The poor agreement between the fit and 
the temperature-dependent data may indicate that the 
two-band model used here still provides an oversimplified 
description of the superconducting gap structure. Thus 
the possible presence of nodes in the gap or of a third 
gap of very low energy cannot be ruled out. 

For a better understanding of the double-gap scenario, 
we consider in detail the frequency dependence of the 
optical conductivity tJi, as evaluated by using the best 
fit parameter values. The model uin and ais of the A- 
and B-band are reported separately in Fig. up to 220 
cm""'^, i.e. in the THz region where the effect of the su- 
perconducting gaps is dominant. In the normal state, 
since ctqa ~ 5500 Q,^^ cmT^ is roughly three times larger 
than (Tqb and 7b is at least ten times larger than 7^1, 
the spectral shape of the total aiN=<yiN A+<yiN b in the 
THz region is dominated by the A-band contribution, as 
shown in Fig. [3]d. Our results are thus in agreement with 



recent findings showing that the normal state conductiv- 
ity of a number of 122 systems [SOj, and in particular 
of BaFe2-£cCo2;As2 compounds (SJl, is given by a sharp 
component superimposed to a much broader component. 
In the superconducting state, since (Tisa (s'isb) vanishes 
below 2A^ (2Ab), the total (Jis=cfisA+'yiSB vanishes 
below the smaller optical gap 2Aa, as shown in Fig. ^p. 
Since, in general, Rs approaches 1 when uis vanishes, 
this explains why only the effect of the smaller gap /S.a 
is well evident in the Rs/Rn spectrum. 

The gap values we obtained are compared in Table 
H] with those reported in recent infrared/THz studies of 
BaFe2-a:Coa:As2 samplcs. The reflectivity measurements 
of van Heumen et al. on a a;=0.14 single crystal jl^ pro- 
vide evidence of a high-energy gap (As =56 cm~^) in very 
good agreement with our finding, while the low-energy 
gap value A^=26.5 cm~^ is higher than our own value 
(15.5 cm~^). Kim et al. [2^ performed reflection and 
ellipsometry measurements on a a:=0.13 single-crystal. 
According to their analysis, three isotropic gaps are nec- 
essary to obtain a good description of the experimental 
data. The low-energy gap Aa value is close to 25 cm~^, 
while the two gaps at higher energy are around 40 and 78 
cm~^. Wu et al. (soj made reflectivity measurements on 
a x=0.16 crystal and their analysis provided evidence of 
a 25 cm~^ gap and also suggested the presence of a sec- 
ond gap at very low energy around 9 cm~^. Gorshunov 
et al. :31], besides reflectance, performed transmittance 
measurements at very low frequencies on a fllm (l7| with 
a;=0.2, and found evidence of only one isotropic gap at 
15 cm~^, in very good agreement with our Aa value. 
THz conductivity spectroscopy measurements by Naka- 
mura et al. [32] on a x=0.2 fllm showed the presence of 
a gap at 22.4 cm"-'^. The reflectivity measurements of 
Nakajima et al. on single crystals gave evidence of 
a single gap at 40 cm~^ for 2;=0.12, around 25 cm~^ for 
x=0.16. Finally, Fisher et al. [s^ measured the complex 
dynamical conductivity of a x=0.2 fllm and their analysis 
gave a nodeless gap at 24 cm~^ and a nodal gap around 
64 cm-i. 

Several different measurements therefore indicate a low 
energy gap Aa mostrly between 15 and 25 cm~^, as re- 
ported in Table ID In trying to explain this uncertainty, 
it might be necessary to consider various effects, such as 
disorder [1, [12, [21], and doping level [81] which can 
differently affect the various samples, and, in the case of 
films, the strain imposed by the substrate T?']. Further- 
more, in the case of standard refiectance measurements, 
it is important to keep in mind that, in terms of signal to 
noise ratio, the quality of THz measurements on a large 
surface film can be better than that of measurements on 
crystals, especially when the wavelength approaches the 
crystal size. The difference among the observed A^ val- 
ues might then be explained, but more bothersome is the 
disagreement in the number and values of the high-energy 
gaps, as summarized in Table [T| 
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TABLE I: THz supeconductings gaps in BaFe2-iCoa; As2 compounds from the present work and literature [2^43^. 



Present work Rcf. [28] Rcf. [29] Ref. [30]" Ref. [31] Ref. [32]'' Rcf. [33] Rcf. [33] Rcfj34] 



Co-doping 


0.16 


0.14 


0.13 


0.16 


0.2 


0.2 


0.12 


0.16 


0.2 


Tc (K) 


22.5 


23 


24.5 


25 


20 


19.9 


25 


20 


25 


N. of gaps 


2 


2 


3 


2 


1 


1 


1 


1 


2 


A A (cm-i) 


15.5 


26.5 


25 


25 


15 


22.4 


40 


25 


24 


Ab (cm-i) 


55 


56 


40 












64 


Ac (cm-l) 






78 















" The authors suggest the presence of a low energy gap around 9 cm ^ , below the range in which reliable data could be acquired. 
^ The authors do not exclude an high energy gap since their measurements are perfomed up to 1 THz (33 cm~^). 



IV. CONCLUSIONS 

We measured the THz reflectance properties of a high 
quality epitaxial film of BaFei.84Coo.i6As2 with Tc=22.5 
K by measuring the R{T)/Rn spectrum, i.e. the reflec- 
tivity ratio between the superconducting state (for 6<T< 
22 K) and the 25 K (normal state) reflectance. The mea- 
sured Rs/Rn spectrum (with Rs = R{6K) first sug- 
gests the presence of a superconducting gap close 
to 15 cm~^. We find that only a two-band model, with 
a gap opening in both bands, can describe the measured 
Rs I Rn spectrum, and in particular its smoothly decreas- 
ing slope. The low-energy gap (Aa=15.5±0.5 cm~^) 
and the high-energy gap Aq (As=55±7 cm^^) corre- 
spond to 2A/kTc values close to 2 and 7, respectively. 
The Rs/Rn data are thus compatible with a two-band 
scenario, with a nodeless isotropic gap opening in both 
bands. Nevertheless, the poor agreement between the 
model and experimental R{T)/Rm at 12 and 18 K can 
be interpreted as an indication for ungapped states at 
low energies. Their presence may be attributed either 
to nodes or to the presence of gaps at energies so low 
as to get thermally excited carriers at relatively modest 
temperatures. 

Despite some differences, the magnitude of the low- 
energy optical gap recently found by different groups sub- 



stantially agree with our A^ value, but there is a large 
disparity on the reported number of gaps and on their 
values. In order to resolve this important issue, THz 
transmission measurements on BaFe2-xCoa;As2 films, to 
be performed in the entire THz region where gap effects 
are expected, could be of great help in settling this de- 
bate. Indeed, in a multigap superconductor, transmit- 
tance measurements are more sensitive to the features re- 
lated to the high-energy gap(s), while reflectance is more 
dependent on the low-energy one, as we have shown in 
the VaSi case [Hi. 

Acknowledgements 

Work at NHMFL was supported under NSF Cooper- 
ative Agreement DlVIR-0084173, by the State of Florida, 
and by AFOSR grant FA9550-06-1-0474. Work at the 
University of Wisconsin was supported by the US De- 
partment of Energy, Division of JVIaterials Science, under 
Award No. DE-FG02-06ER463. The authors thank M.S. 
Rzchowski for helpful discussion. Work at the University 
of Rome was partially funded by CARIPLO Foundation 
(Project 2009-2540 "Chemical Control and Doping Ef- 
fects in Pnictide High-temperature Superconductors"), 
at the University of Genova by the "Italian Foreign Af- 
fairs Ministry (IMAE) - General Direction for the Cultural 
Promotion" . 



[1] Kamihara Y., Watanabe T., Hirano M. and Hosono H., 
J. Am. Chem. Soc. 130, (2008) 3296 

[2] Qazilbash M.M., Hamlin J. J., Baumbach R.E., Lijun 
Zhang, Singh D.J. , Maple M.B. and Basov D.N., Na- 
ture Physics 5, (2009) 647 

[3] Yang W.L., Sorini A.R, Chen C-C, Moritz B., Lee W.- 
S., Vernay F., Olalde-Velasco P., Denlinger J.D., Delley 
B., Chu J.-H., Analytis J.G., Fisher LR., Ren Z.A., Yang 
J., Lu W., Zhao Z.X., van den Brink J., Hussain Z. , 
Shen Z.-X. and Devereaux T.P., Phys. Rev. B 80, (2009) 
014508 

[4] Muschler B., Prestel W., Hackl R., Devereaux T.P., An- 
alytis J.G., Jiun-Haw Chu and Fisher I.R., Phys. Rev. B 
80, (2009) 180510(R) 



[5] Goko T., Aczel A.A.,Baggio-Saitovitch E., Bud'ko S.L., 
Canfield P.C., Carlo J.P., Chen G.F.,Dai P., Hamann 
A.C.,Hu W.Z., Kageyama H., Luke G.M., Luo J.L., 
Nachumi B., Ni N., Reznik D., Sanchez-Candela D.R., 
Savici A.T., Sikes K.J., Wang N.L., Wiebe C.R., 
Williams T.J., Yamamoto T., Yu W. and Uemura Y.J., 
Phys. Rev. B 80, (2009) 024508 

[6] Salem-Sugui Jr. S. , Ghivelder L., Alvarenga A.D., Pi- 
mentel Fr. J.L., Luo H., Wang Z. and Wen H.H., Phys. 
Rev. B 80, (2009) 014518 

[7] Gordon R.T., Ni N., Martin C., Tanatar M.A., Vannette 
M.D., Kim H., Samolyuk G.D., Schmalian J., Nandi S., 
Kreyssig A., Goldman A. I., Yan J.Q., Bud'ko S.L., Can- 
field P.C., Prozorov R., Phys. Rev. Lett. 102, (2009) 



6 



127004 

[8] Terashima K., Sekiba Y., Bowen J. H., Nakayama K., 
Kawahara T., Sato T., Richard P., Xu Y.-M., Li L. J., 
Cao G. H., Xug Z.-A., Ding H. and Takahashi T., Proc. 
Natl. Acad. Sci. 106, (2009) 7330 
[9] Samuely P., Pribulova Z., Szabo P., Pristas G., Budko 
S.L. and Canfield P.C., Physica C 469, (2009) 507 

[10] Williams T.J., Aczel A. A., Baggio-Saitovitch E., Bud'ko 
S.L., Canfield P.C., Carlo J. P., Goko T., Munevar J.,Ni 
N., Uemura Y.J., Yu W. and Luke G.M., Phys. Rev. B 
80, (2009) 094501 

[11] Hardy P., Wolf T., Fisher R.A., Eder R., Schweiss P., 
Adelmann P. , v. Lohneysen H. and Meingast C, Phys. 
Rev. B 81, (2010) 0060501 

[12] Tanatar M.A., Reid J. P., Shakeripour H., Luo X.G., 
Doiron-Leyraud N., Ni N., Bud'ko S.L., Canfield P.C., 
Prozorov R. and Taillefer L., Phys. Rev. Lett. 104, 
(2010) 067002 

[13] Luan L., Auslaender O.M., Lippman T.M., Hicks C.W., 
Kalisjy B., Chu J.H., Alalytis J.G., Fisher LR., Kirtley 
J.R. and Moler K.A., Phys. Rev. B 81, (2010) 100501 

[14] Mazin I.L, Singh D.J., Johannes M.D. and Du M.H., 
Phys. Rev. Lett. 101, (2008) 057003 

[15] Luo H., Wang Z., Yang H., Cheng P., Zhu X. and Wen 
H.H. Supercond. Sci. Technol. 21, (2008) 125014 

[16] Lee S., Jiang J., Nelson C.T., Bark C.W., Weiss J.D., 
Tarantini C, Jang H.W., Folkman CM., Back S.H., 
Polyanskii A., Abraimov D., Yamamoto A., Zhang Y., 
Pan X.Q., Hellstrom E.E., Larbalestier D.C. and Eom 
C.B. Nature Materials 9, (2010) 397 

[17] lida K., Hanisch J., Hiinhe R., Kurth F., Kidszun M., 
Haindl S., Werner J., Schultz L. and Holzapfel B., Appl. 
Phys. Lett. 95, (2009) 192501 

[18] Sefat A.S., Jin R., McGuire M.A., Sales B.C., Singh D.J. 
and Mandrus D., Phys. Rev. Lett. 101, (2008) 117004 

[19] Yin Y., Zech M., Williams T.L., Wang X.F., Wu G., 
Chen X.H. and Hoffman J.E. Phys. Rev. Lett. 102, 

(2009) 097002 

[20] Massee F., Huang Y., Huisman R., de Jong S., Goed- 
koop J.B. and Golden M.S., Phys. Rev. B 79, (2009) 
220517(R) 

[21] Hashimoto K., Shibauchi T., Kasahara S., Ikada K., 
Tonegawa S., Kato T., Okazaki R., van der Beek C.J., 
Konczykowski M., Takeya H., Hirata K., Terashima T. 
and Matsuda Y., Phys. Rev. Lett. 102, (2009) 207001 

[22] Carbotte J. P. and Schachinger E. Phys. Rev. B 81, 

(2010) 104510 

[23] Burns G. Solid state physics (Academic Press, Boston, 
1990) 



[24] Mattis D.C. and Bardeen J., Phys. Rev. Ill, (1958) 412 
[25] Tinkham M. Introduction to Superconductivity 

(McGraw-Hill, New York, 1975) 
[26] Hu W.Z., Zhang Q.M. and Wang N.L. Physica C 469, 

(2009) 545 

[27] Drechsler S.L., Rosner H., Grobosch M., Behr G., Roth 
F., Fuchs G., Koepernik K., Schuster R., Malek J., El- 
gazzar S., Rotter M., Johrendt D., Klauss H.H., Biichner 
B. and Knupfer M. arXiv:0904.0827 preprint, 2009 

[28] van Heumen E., Huang Y., de Jong S., Kuzmenko A.B., 
Golden M.S. and van der Marel D., arXiv:0912.0636l 
preprint, 2009 

[29] Kim K.W., Rossle M., Dubroka A., Malik V.K., Woff T. 
and C. Bernhard, arXiv:0912.0140 preprint, 2009 

[30] Wu D., Barisic N., Kallina P., Faridian A., Gorshunov B., 
Drichko N., Li L.J., Lin X., Cao G.H., Xu Z.A., Wang 
N.L. and Dressel M., Phys. Rev. B 81, (2010) 100512(R) 

[31] Gorshunov B., Wu D., Voronkov A. A., Kallina P., lida 
K., Haindl S., Kurth F., Schuhz L., Holzapfel B. and 
Dressel M., Phys. Rev. B, 81 (2010) 060509(R) 

[32] Nakamura D., Imai Y., Maeda A., Katase T., Hiramatsu 
H., and Hosono H., larXiv: 0912.4351v2 preprint, 2010 

[33] Nakajima M., Ishida S. , Kihou K., Tomioka Y., Ito T., 
Yoshida Y., Lee C. H., Kito H., lyo A., Eisaki H., Kojima 
K. M., and Uchida S., Phys. Rev. B 81, (2010) 104528 

[34] Fischer T., Pronin A. V., Wosnitza J., lida K. , Kurth 
F. , Haindl S., Schultz L., Holzapfel B., and Schachinger 

E. , arXiv;1005.0692vl preprint, 2010. 

[35] Lupi S., Nucara A., Perucchi A., Calvani P., Ortolani 
M., Quaroni L. and Kiskinova M., J. Opt. Soc. Am. B 
24, (2007) 959 

[36] Ortolani M., Dore P., Di Castro D., Perucchi A., Lupi S., 
Ferrando V., Putti M., Pallecchi L, Ferdeghini C. and Xi 
X.X., Phys. Rev. B 77, (2008) 100507(R) 

[37] Kuzmenko A.B.. Phys. C 456, (2007) 63 

[38] Perucchi A., Nicoletti D., Ortolani M., Marini C, Sopra- 
case R., Lupi S., Schade U., Putti M., Pallecchi L, Taran- 
tini C, Ferretti M., Ferdeghini C, Monni M., Bernardini 

F. , Massidda S. and Dore P., Phys. Rev. B 81, (2010) 
092509) 

[39] Zimmermann W., Brandt E.H., Bauer M., Seider E. and 

Genzel L., Phys. C 183, (1991) 99 
[40] Dressel M. and Griiner G., Electrodynamics of Solids 

(Cambridge University Press, Cambridge, England, 

2002) 

[41] Dore P., Paolone A., and Trippetti R., J. Appl. Phys. 80, 
(1996) 5270; Dore P., De Marzi G. and Paolone A., Int. 
J. of ir and mm Waves 18, (1997) 125 



